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Abstract 

We investigate neutrino-nucleus interactions at energies around 1 GeV. In this 
regime, the main contributions to the cross sections come from quasi-elastic and 
A production processes. Our formalism, based on the Impulse Approximation is 
well suited to describe both types of interactions. We focus on a series of important 
nuclear effects in the interaction of electron neutrinos with 16 O, also relevant to 
future /3-Beams. Our results show that the Fermi gas model, widely used in data 
analysis of neutrino experiments, overestimates the total cross sections by as much 
as 20 %. 
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1 Introduction 



The field of neutrino physics is rapidly developing after atmospheric and solar 
neutrino oscillations have been established. The results of atmospheric, solar, 
accelerator and reactor neutrino experiments |T] show that flavour mixing 
occurs not only in the hadronic sector, as it has been long known, but in the 
leptonic sector as well. The experimental results point to two very distinct 
mass difference^], Am^ ps 8.2 x 1CT 5 eV 2 and |Am 2 tm | fa 2.5 x 10~ 3 eV 2 . 
Only two out of the four parameters of the three-family leptonic mixing matrix 
Upmns [I] are known: 9 12 ~ 34° and #23 ~ 45°. The other two parameters, #13 
and S, are still unknown: for the mixing angle #13 direct searches at reactors [5] 
and three-family global analysis of the experimental data [6fT] give the upper 

1 A third mass difference, /±m 2 LSND ~ 1 eV 2 , suggested by the LSND experiment 
[2], has not been confirmed yet [3]. 
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bound #13 < 11.5°, whereas for the leptonic CP-violating phase S we have no 
information whatsoever. Two additional discrete unknowns are the sign of the 
atmospheric mass difference and the #23-octant (if # 2 3 7^ 45°). 

The full understanding of the leptonic mixing matrix, together with the dis- 
crimination of the Dirac/Majorana character of neutrinos and with the mea- 
surement of their absolute mass scale, represents the main goal of neutrino 
physics in the next decade. The SK detector has gathered indirect evidence 
of — > v T conversion of atmospheric neutrinos whereas the SNO detector 
[H] has shown that a fraction of the z/ e 's emitted by the Sun core reaches the 
Earth-located detectors converted into v^s and z/ T 's (and not into unobserv- 
able sterile neutrinos). On the other hand, new-generation experiments have 
been proposed to look for the fleeting and intimately related parameters #13 
and 5 through the more promising "appearance channels" such as v e — > (the 
"golden channel" [S]) and v e — > v T (the "silver channel" [ID])- However, strong 
correlations between #13 and 5 and the presence of parametric degeneracies in 
the (613, S) parameter space [TT][T2"] . make the simultaneous measurement of 
the two variables extremely difficult. A further problem arises from our present 
imprecise knowledge of atmospheric parameters, whose uncertainties are far 
too large to be neglected when looking for such tiny signals as those expected 
in appearance experiments [13]. Most of the proposed solutions to avoid all 
these problems suggest to combine different experiments and facilities, such 
as Super-Beams (of which T2K [H] is the first approved one), /3-Beams [15] 
or Neutrino Factories [IB], although it is not completely clear whether the 
sensitivity of such experiments will be enough. 

In view of these developments, it is vital to reduce as much as possible any 
source of uncertainty that could ruin the required precision and, among them, 
the knowledge of the neutrino-nucleus cross sections is one of the thorniest 
problems to be faced. In fact, the data are either very few (the case of neutri- 
nos) or missing altogether (the case of antineutrinos). On top of that, the few 
available data have generally not been taken on the same targets used in the 
experiments (either water, iron, lead or plastics), and the extrapolation from 
different nuclei is complicated by nuclear effects, that can play an important 
role at the considered energies. Among the proposed experiments, Minerva 
[T7] and SciBooNE [THj will have the capability for precise measurements in a 
wide range of energies and different nuclear targets. 

In this paper we will focus on total nuclear cross sections at energies up to ~ 
1 GeV, relevant to a number of neutrino experiments as well as to the planned 
Super-Beams and /5-Beams. As an example of the latter, in Fig. [U we show 
the energy spectra for two different setups, in which the three decay modes of 
18 Ne [21] are considered. 

In this energy regime, the dominant contribution to the charged lepton produc- 
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Fig. 1. (3-Beam fluxes as a function of neutrino energy for two specific sets 
of 7 values. Low and high 7 correspond to (7^,7^) = (100,60) /TP) /, and 
(7f.,7i0 = (250,150) JH, respectively. 

tion cross section comes from quasi-elastic reactions. However, for E v ^ 0.5 GeV, 
inelastic production of charged leptons through excitation of the A (as well 
as other) resonance becomes relevant. At higher energies, deep inelastic scat- 
tering also becomes progressively important, but the impact on the energy 
regime under investigation is expected to be negligible. Therefore, we do not 
take this contribution into account. 

Quasi-elastic scattering, pion production and deep inelastic scattering are com- 
bined in event generators (e.g. NUANCE [22], NEUGEN [23] and NEUT [23]). 
that model the neutrino detector response using a simplified picture, referred 
to as Relativistic Fermi Gas Model (RFGM), in which the nucleus is described 
as a collections of quasi-free nucleons. A wealth of high-precision electron scat- 
tering data have shown that the RFGM neglects important features of nuclear 
dynamics [25]. The most important is the presence of strong nucleon-nucleon 
correlations, leading to the appearance of high momentum and high energy 
components in the energy-momentum distribution of the nucleons in the target 
nucleus [2T>] . 

In this paper we go beyond the simple RFGM, and analyse neutrino-nucleus 
interactions at intermediate energies using a more realistic description of nu- 
clear dynamics, based on nuclear many-body theory (NMBT). We evaluate 
the cross sections for inclusive processes, in which only the outgoing charged 
lepton is detected, including both quasi-elastic and pion production processes. 
Our formalism is based on the impulse approximation (IA) picture, which 
is expected to be applicable at the large momentum transfers corresponding 
to beam energies around 1 GeV. Final state interactions (FSI) between the 
hadrons produced at the elementary weak interaction vertex and the specta- 
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tor nucleons have been not taken into account, as they do not contribute to 
the total cross section. However, the effect of statistical correlations between 
the struck nucleon and the spectators, which are known to be important in 
quasi-elastic scattering at low Q 2 , has been included. 

The plan of the paper is the following. In Section [2] we recall the formalism 
of the IA used to describe charged current (CC) interactions. In Section [3] 
we show our results for the inclusive CC cross sections, while in Section H] we 
present the results for resonance production. Finally we draw the conclusions 
and outline the future prospects of our work. The appendices contain formulae 
that, together with those presented in the main text, will enable the readers 
to write their own codes to compute the cross sections discussed in this paper. 



2 The Impulse Approximation 



The differential cross section for the process + A — > t~ + X (Fig. [2D in which 
a neutrino carrying initial four- momentum k = (E u , k) scatters off a nuclear 
target to a state of four- momentum k = (E#, k ), the target final state being 
undetected, can be written in Born approximation: 



d 2 a G|K 2 d |k, 



dVtdE, 16 7T 2 Ikl 



V^T, (i) 



where Gp is the Fermi constant and V u d is the CKM matrix element coupling 
u and d quarks. The leptonic tensor, that can be written in the form 



k„ k v + k v k^ (j^y{k ■ k 



IE 



k'P k a 



(2) 



is completely determined by lepton kinematics, whereas the nuclear tensor 
W^ u , containing all the information on strong interactions dynamics, describes 
the response of the target nucleus. Its definition involves the initial and final 
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Fig. 2. Feynman diagram for the process vg + A 
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hadronic states |0) and \X), carrying four momenta po and px, respectively, 
as well as the nuclear electroweak current operator J\: 



wr 



£<0|J^|X) (X\r A \0) SMfa + q-px) 



(3) 



where the sum includes all hadronic final states. Calculation of W^ u at mod- 
erate momentum trasfers (|q| <0.5 GeV) can be carried out within NMBT 
using nonrelativistic wave functions to describe the initial and final states and 
expanding the current operator in powers of q/m^v, rn^ being the nucleon 
mass. However, at higher values of |q|, corresponding to E v > 1 GeV, we can 
no longer describe the final states \X) in terms of nonrelativistic nucleons. 
Calculation of i n this regime requires a set of simplifying assumptions, 
allowing one to take into account the relativistic motion of final state parti- 
cles earring momentum ~ q, as well as the occurrence of inelastic processes, 
leading to the appearance of hadrons other than protons and neutrons. 



In the rest of the paper, we adopt the IA scheme, based on the assumptions 
that at large enough q the target nucleus is seen by the probe as a collection of 
individual nucleons and that the particles produced at the interaction vertex 
and the recoiling (A — l)-nucleon system evolve independently (see Fig. [3]for a 
pictorial representation of the IA picture). As a consequence, IA neglects both 
statistical correlations due to Pauli blocking and the rescattering processes 
driven by strong interactions (FSI). 




Fig. 3. Diagrammatic representation of the process V£ + A 
Approximation scheme. 



+ X in the Impulse 



Within this picture, the nuclear current can be written as a sum of one-body 
currents, i.e. J A — > J2i Ji, while the final state reduces to the direct product 
of the hadronic state produced at the weak vertex (with momentum p ) and 
that describing the (A — l)-nucleon residual system, with momentum p^: 
\X) — > |i,p) <8> \R., Pii)- The calculation of the weak tensor is described in 
Ref. [27]. Here we only quote the final result: 



^A=lJd 3 pdEP(p,E) 



4E\ P \ £| p+q | 



W^(p,q) 



(4) 
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where E p = J |p| 2 + m 2 N . The above definition deserves some comments. The 
function P(p, E) is the target spectral function, i.e. the probability distribution 
of finding a nucleon with momentum p and removal energy E in the target 
nucleus. It then encodes all the informations about the initial (struck) particle. 
The quantity W^ u is the tensor describing the weak interactions of the z-th 
nucleon in free space; the effect of nuclear binding of the struck nucleon is 
accounted for by the replacement: 

q = 0> q) ->• q = iy, q) 

with v = E\ p+q \ — E\ p \. In particular, the variable v takes into account the 
fact that a fraction 5v of the energy transfer goes into excitation energy of the 
spectator system, and that the elementary scattering process can be described 
as if it took place in free space with energy transfer v = v — 8v. It follows 
that the second argument in the hadronic tensor is p = (E\ p \, p). Substituting 
Eq. PJ) into Eq. (JTJ, we get the final formula for the nuclear cross section: 



d2<TlA f d 3 P dEP(p,E)^^, (5) 



dndE e J r y ^ ' dfldE L 
in which we have redefined the elementary cross section as 

d 2 0~elem _ Gp V^ d |k | 1 ^ . . 

cKldEt 32 vr 2 |k| AE p E\ p+fl \ 1 ' 

The integration limits for Eq. (jHJ) are given in Appendix lAl 

The hadronic tensor W^ v can be written in terms of five structure functions 
Wi as 



Wo a W-i Wa 

W» v = -<T Wx+FF^l + i q a p p -f- + q» <f + (7) 



77% m N m N 



Wk 



m 



N 



Its contraction with the leptonic tensor (j2J) yields 



IT W, v = 16 £ Wi ( 4-1 • (8) 
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N , 



The functions Ai, containing kinematical factors, are collected in Appendix 
IB1 whereas the structure functions Wi, expressed in terms of form factors, are 
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given in Appendix O 

Eq. (J3J) is a general formula that can be used to describe quasi-elastic scattering 
and resonance production in neutrino-nucleus charged current interactions. 
The explicit expressions differ in the analytical form of the relevant form 
factors and for the replacement of the energy conserving 5-function with a 
Breit-Wigner factor accounting for the finite width of the resonance. 



3 Elastic cross sections 

To obtain the elastic cross sections from Eq. we need to clarify the form of 
the structure functions Wi and the meaning of the spectral function. The first 
part of the job is easily done once we express the structure functions in terms 
of form factors F\, F2, Fa, Fp using the hypotheses of CVC (Conserved Vector 
Current) and PCAC (Partially Conserved Axial Current), as usually done 
for scattering off a free nucleon. In addition, the definition of the structure 
functions includes an energy conserving 5-function 5(s — m 2 N ), where s = 
(p + q) 2 is the squared invariant mass of the final state of the struck nucleon. 
In Appendix O we collect the expression of Wi in terms of form factors. Note 
that, except where differently stated, we use the values My = 0.84 GeV and 
M A = 1.05 GeV for the vector and axial mass, respectively. 



3. 1 Fermi gas and spectral function 

The RFGM [28] , widely used in Monte Carlo simulations, provides the simplest 
form of the spectral function: 

P RFGM (p, E) = \ —^ \ 0(Pf ~ p) S(E p — E B + E) , (9) 

where pp is the Fermi momentum and E B is the average binding energy, in- 
troduced to account for nuclear binding. The term in parenthesis is a constant 
needed to normalize the spectral function to the number of target nucleons, 
A. Thus, in this model pp and E B are two parameters that are adjusted to 
reproduce the experimental data. For oxygen, the analysis of electron scat- 
tering data yields pp = 225 MeV and Ep = 25 MeV [29]. In more refined 
versions of RFGM, the binding energy of the nucleon is not a constant but a 
momentum-dependent function V(|p|) [30], entering the 5-function in Eq. (Q 
according to 5(E P + V(\p\) + E). 
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Electron scattering data have provided overwhelming evidence that the energy- 
momentum distribution of nucleons in the nucleus is quite different from that 
predicted by RFGM (see, e.g. Ref. [31]). The most important feature emerged 
from these data is the presence of strong nucleon-nucleon (NN) correlations. 

Strong dynamical correlations give rise to virtual scattering processes leading 
to the excitation of the participating nucleons to states of energy larger than 
the Fermi energy, thus depleting the single particle levels within the Fermi sea. 
As a consequence, the spectral function associated with nucleons belonging to 
correlated pairs extends to the region of |p| 3> Pf and E S> E B . 

Clearcut evidence of correlation effects has been provided by electron- and 
hadron-induced nucleon knock-out experiments [52]. While the spectroscopic 
lines corresponding to knock out from the shell model states are clearly seen in 
the measured energy spectra, the associated strengths are consistently lower 
that expected, regardless of the nuclear mass number. Unfortunately, system- 
atic measurements of the missing strength (typically about 20%) pushed to 
higher energy by NN correlations, have not been carried out yet. However, the 
results of a pioneering JLab experiment [33] appear to be consistent with the 
expectations based on lower energy data. 

Correlation effects can be consistently taken into account within NMBT, in 
which nuclei are described in terms of nonrelativistic nucleons interacting 
through the hamiltonian 

H a = £^+ E % + ..., (10) 
i=l z " i j>i=l 

where Vij provides a quantitative account of the properties of the two-nucleon 
system, i.e. deuteron properties and ~ 4000 accurately measured nucleon- 
nucleon scattering phase shifts at energies up to the pion production thresh- 
old [31], while the ellipses refer to a small additional three- nucleon potential. 
The energies of the ground and low-lying excited states of nuclei with mass 
number A < 10 calculated within NMBT are in excellent agreement with the 
experimental values |35j . 

The calculation of P(p, E) within NMBT involves a degree of complexity 
that rapidly increases with A, so that it has been only carried out for A < 
4 [37)f3"T|r3 8 39 4 0f4"T] . However, thanks to the simplifications associated with 
translation invariance, highly accurate results are also available for uniform 
nuclear matter, i.e. in the limit A — > oo with Z=A/2 [4"2l43j {Z denotes the 
number of protons). 

The spectral functions for medium-heavy nuclei, ranging from Carbon {A = 
12) to Gold (A = 197), have been modeled using the Local Density Approx- 
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imation (LDA) pH] , i n which the experimental information obtained from 
nucleon knock-out measurements is combined with the results of theoretical 
calculations of the nuclear matter P(p, E) at different densities. Within this 
approach, P(p, E) is divided in two parts, corresponding to low momentum 
nucleons, occupying shell model states, and high momentum nucleons, respec- 
tively. 

The correlation contribution to P(p, E) of uniform nuclear matter has been 
calculated by Benhar et al. for a wide range of density values [H]. Within the 
LDA scheme, the results of Ref. [H] can be used to obtain the corresponding 
quantity for a finite nucleus of mass number A. The full LDA nuclear spectral 
function can then be written as 



From the spectral function, we can obtain the nucleon momentum distribution, 
defined as 



The nucleon momentum distributions of 16 O and 197 Au obtained from LDA, 
normalized to unity, are shown in Fig. HI For reference, the RFGM momentum 
distribution corresponding to Fermi momentum pp = 225 MeV is also shown 
by the dashed line. 

The results of Fig. [4] show that n{p)/A becomes nearly independent of A at 
large |p| 300 MeV/c). This feature suggests that the correlation part of 
the spectral function also scale with the target mass number, so that nuclear 
matter results can be used at finite A. It also clearly appears that the LDA 
n(p) is very different from that obtained from the RFGM, in which the large 
momentum tail is missing. 

As pointed out in Section [TJ within IA all FSI are neglected. The validity 
of this approximation has been thoroughly analyzed within the framework of 
non relativistic models of the nuclear response using correlated basis function 
(CBF) perturbation theory [12]. 

In quasi-elastic inclusive processes dynamical FSI, which are long known to 
be important in semi-inclusive and exclusive electron-nucleus scattering pro- 
cesses, are much weaker. Their main effects are i) an energy shift of the differ- 
ential cross section, due to the fact that the struck nucleon feels the mean field 
generated by the spectator particles and ii) a redistribution of the strength, 
leading to the quenching of the peak and the enhancement of the tails. 

Within non relativistic NMBT, it can be shown that FSI do not affect the total 



P LDA (P, E) = P MF (p, E) + P corr (p, E) . 



(11) 




(12) 
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Fig. 4. Momentum distribution ofnucleons in the oxygen (solid line) and gold (dotted 
line) ground states. Dashed line: RFGM with Fermi momentum pf = 225 MeV and 
E B = 25 Me V. 

inclusive cross section, resulting from integration over the lepton variables 
(energy loss and scattering angle). This result follows from the fact that the 
final states of the A-body system can be generated from the complete set of 
asymptotic states, corresponding to independent propagation of the struck 
nucleon and the recoiling spectator system, through a unitary transformation. 
When the energy of the knocked out nucleon is large and relativistic effects 
cannot be neglected, this procedure can still be applied using a generalization 
of Glauber theory, based on the eikonal approximation, as discussed in, e.g., 
Ref. [45J. As our work is focused on total cross sections only, FSI effects will 
not be taken into account. 

The effect of statistical correlations, arising from the antisymmetrization of 
the struck nucleon with respect to the spectator particles, is also disregarded 
within IA. The results of calculations carried out within nonrelativistic models 
[12] show that Pauli blocking has negligible effect on the differential cross 
section for three- momentum transfer |q| larger than ~ 500 MeV. However, as 
in the kinematical regime discussed in our work the total cross section includes 
significant contributions from scattering at lower |q|, the modification of the 
phase space available to the knocked-out nucleon must be taken into account. 
A rather crude prescription to estimate the effect of Pauli blocking is based 
on the replacement [27] 




(13) 
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where p F is the average nuclear Fermi momentum, defined as 
Vf= \ d3r Pa(t)pf(t), 



(14) 



with pp (r) = (37i 2 pA(r)/2) 1 / 3 , pa(t) being the nuclear density distribution. For 
oxygen, Eq. ffT4"|) yields p F = 209 MeV. Note that, unlike the spectral function, 
the quantity defined in Eq. (TT~3T> does not describe intrinsic properties of the 
target only, as it depends explicitely on the momentum transfer. Overall, the 
approach described above involves no adjustable parameters (as compared 
with other models discussed in Section ET21 . 

The relative weight of different nuclear effects is illustrated in Fig. [51 showing 
the total cross section of the process v e + 16 O — > e~ + X as a function of 
neutrino energy. The dot-dash line represents eight times the elementary cross 

quasi-elastic inclusive cross section 

10 




E v (GeV) 

Fig. 5. Total quasi-elastic cross section for the process v e 16 — ► e~ X. The dot-dash 
line represents eight times the elementary cross section; the dashed line is the result 
of the RFGM with Fermi momentum pp = 225 MeV and binding energy Eb = 25 
Me V; the dotted and solid lines have been obtained using the LDA spectral function, 
with and without inclusion of Pauli blocking, respectively. 

section, i.e. the cross section for interactions of neutrinos on free neutrons. The 
dashed line is the result of the RFGM without inclusion of Pauli blocking (PB). 
The difference between the dot-dash and dashed lines, the latter being sizably 
smaller, results from the Fermi motion (with momenta not exceeding pp) of 
the target nucleons. The solid line corresponds to our results obtained from 
the LDA oxygen spectral function. It clearly appears that the inclusion of NN 
correlations leads to a further reduction of the cross section. This reduction 
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is further enhanced if we include PB, as shown by the dotted line, since the 
available phase space for the knocked nucleon gets smaller, see Eq. ( !T3|) . Note 
that, as expected, the importance of PB decreases as neutrino energy increases. 
For example, the difference between the solid and dotted lines turns out to be 
10 % and 6 % at neutrino energy 0.5 GeV and 1.5 GeV, respectively. It has to 
be mentioned that the prescription of Eq. (fT3l) can also be used to include PB 
in the RFGM. The size of the effect is the same as in the calculation using the 
spectral function, the resulting cross section being suppressed by ~ 5 — 10 % 
with respect to the dashed line of Fig. [5J Overall, based on the results shown 
in Fig. [5] one can safely state that a realistic treatment of neutrino-nucleus 
interactions must go beyond the RFMG. 

The same conclusions can be drawn for antineutrino interactions. Neglecting 
the isospin breaking effect leading to the proton to neutron mass difference, 
the same Eqs. (jSJ)- ([6]) can be applied. The only difference arises in the sign in 
front of the parity- violating term of Eq. ([7J), £ MJ , Qj g, which should be replaced 
by —s^vap- Therefore, we do not show the corresponding figure. 



3.2 Comparison with other models 

Over the past few years, many authors have investigated neutrino cross sec- 
tions and the related nuclear effects (see, e.g., Ref. [IS]). In this section we 
compare our results with the existing literature for the case of oxygen and 
carbon targets. It has to be kept in mind, however, that, while many of the 
available results correspond to E v < 100 MeV [IT], the IA scheme adopted 
in our work is not expected to be accurate in this regime. Therefore, our 
comparison is mostly focused on the region of higher energies. 

In Ref. [IB] , the CC quasi elastic cross section is obtained from the self-energy 
of the gauge bosons W in the nuclear medium, whose calculation is carried 
out using the medium-modified nucleon propagator resulting from a LDA 
implementation of the Fermi gas model. The effects of long range nuclear 
correlation (RPA) are also included by means of an effective NN interaction 
of the Landau-Migdal type. Coulomb distortion effects, accounting for the 
fact that the charged lepton produced in neutrino interaction is moving in 
the Coulomb field of the nucleus, are also taken into account. A collection 
of the resulting <T( 16 0(V e , e~)X) is reported in Tab. IV of Ref. [H], in which 
relativistic and non-relativistic nucleon kinematics, as well as the inclusion of 
FSI and RPA effects are analyzed. 

In Fig. [6] we compare our results to those of Ref. [H]. The solid line corresponds 
to our calculation with the LDA spectral function and PB included; the shaded 
area is bounded by the minimum and maximum values of the cross section 
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of Ref . [48j . For reference, the RFGM prediction is also included. It can be 
seen that the agreement between our results and those reported in Ref. |48] is 
satisfactory. The small discrepancy at low neutrino energy, below ~ 0.15 GeV, 
is likely to be ascribed to the inadequacy of the IA to describe small momentum 
transfer processes. It is apparent thet the RFGM largely overestimates the 
cross section. 

At higher energies, we compare our results with those of Refs. [49j and |50j . 
In the first paper, <r( 16 0(y e , e~)X) is calculated using the same formalism as 
described in [48]. The corresponding results are presented in the top panel of 

Fig.m 

In order to pin down the source of the sizable disagreement between our results 
and those of Ref. [49j . we have tested the sensitivity of the cross section to 
the vector and axial masses using two different sets of values: My = 0.86 and 
M A = 1.05, as employed by Ahmad et al. [U], and My = 0.84 and M A = 1.03. 
It clearly appears that the discrepancies at E u > 0.5 GeV cannot be removed 
by a different choice of My and Ma, and should rather be ascribed to the 
different descriptions of nuclear dynamics. 

A much better agreement is observed in the bottom panel of Fig. [TJ in which we 
compare our results for <t( 16 0(^ m , fi~)X), to those obtained in Ref. [50] using 
different models. The approach of Maieron et al, based on a IA scheme similar 
to the one adopted in our work, uses Relativistic Shell-Model wave functions 
to describe bound nucleon states and two different prescriptions to include 
FSI effects. The final state of the knocked out nucleon is described using i) 
distorted waves obtained solving the Dirac equation with a phenomenological 
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Fig. 6. Comparison between our results for o~( ie O(v e , e~)X) (solid line) and those 
presented in J^g| /; the shaded region represents the range of cross sections predicted in 
Ref. 148^ (see text for details). The RFGM results corresponding to Fermi momentum 
Pf = 225 MeV and Eb = 25 MeV are also shown. 
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Fig. 7. Top: comparison between o^O(y e , e~)X) resulting from our calcula- 
tions (solid and dashed lines) and that of Ref. \^9^ . Bottom: comparison of 
rj( 16 0(f^, n~)X) resulting from our calculations and those obtained in Ref. \50\/ us- 
ing different models (see text for details). 



relativistic complex optical potential (ROP) or ii) the continuum solutions of 
the same Dirac equation used for the initial bound states (Relativistic Mean 
Field (RMF) approach. 

The diamonds in the bottom panel of Fig. [7] show the RMF results of Ref. 
[5U] , whereas triangles- up and circles have been obtained by the same authors 
including the real part only and the full ROP (real + imaginary part), re- 
spectively. The most stricking feature is the large effect of the imaginary part 
of the ROP, leading to a suppression of the total cross section of more than 
40 %. This effect can be easily understood, as the imaginary optical potential 
produces a loss of flux in the one particle-one hole channel. However, it con- 
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tradicts the well known result stating that, in the kinematical regime in which 
incoherent scattering dominates, the integrated inclusive cross section is not 
affected by FSI [ST] - In fact, it has long been realized that, unless the reap- 
pearence of the missing flux in the two particle-two hole channel is properly 
taken into account, the prescription employed in Ref. [50J leads to a violation 
of the non energy weighted sum rule of the nuclear response [51] . On the other 
hand, Fig. [7] shows that the cross sections obtained by Maieron et al. neglect- 
ing FSI effects or using only the real part of the ROP, whose effect is very 
small, are in excellent agreement with our results. 



4 Inelastic cross sections 

One pion production is the second process which contributes to the total 
cross section in the kinematical regime discussed in this paper; it proceeds 
mainly through resonance production, the leading contribution coming from 
excitation of the A, i.e. through the the processes 

v e + n^£ + A + , v + p-+£ + A ++ . (15) 

In this section, we calculate the inclusive resonance production cross section 
applying the same formalism used for quasi-elastic interactions, in which nu- 
clear effects in the initial state are described by the spectral function. Obvi- 
ously, PB has no impact on the A produced at the weak interaction vertex, 
although it may affect the phase space available to the nucleon appearing in 
the aftermath of its decay. It has to kept in mind, however, that the inclu- 
sive cross section is only sensitive to processes taking place within a distance 
~ l/|q| of the primary vertex, so that at large enough |q| it is unlikely to be 
affected by two-step reaction mechanisms. 

4-1 A resonance production 

To account for A production, we can apply Eqs. fl5J)-(jSD with only minor 
changes. Unlike the quasi-elastic case, the structure functions depend now on 
two variables, namely q 2 and W 2 , where W 2 is the squared hadronic invariant 
mass. The energy conserving (^-function is replaced by the Breit-Wigner factor 

M R T R 1 

W - M 2 R ) 2 + M 2 R T 2 R ' (16) 

where M R is the resonance mass and T R its decay width, and an additional 
integration in the variable W must be carried out. For the dependence of the 
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structure functions on form factors for A ++ production, we closely follow Ref. 
|52p1 and use isospin symmetry to relate A ++ and A + form factors through 

<A++iJ>) = V3(A + |J>) . (17) 



It has to be mentioned that, in principle, the form factors entering the calcu- 
lation of the cross section could be extracted from proton and deuteron data. 
However, electron scattering studies show that extrapolations of the existing 
phenomenological models to the region of low Q 2 leads to sizably underesti- 
mate the nuclear cross sections in the A production region [271153] . 



Our results for CC A production via 16 0(z/ e , e~) interactions are shown in 
Fig. [HI and compared with the elementary cross section and the RFGM pre- 
diction. 
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A production 
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our result 
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E v (GeV) 
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Fig. 8. CC A production via 16 0(z^ e ,e - ) interactions. The solid line shows the re- 
sults obtained using the LDA spectral function; the dashed and dot-dash curves are 
the elementary cross section and the RFGM prediction (with pp = 225 MeV and 
Eb = 25 MeV), respectively. 

As in the elastic case, nuclear effects strongly reduce the cross section with 
respect to the free space value; moreover, the RFGM appears to overestimate 
A production, an effect that has been already observed by different authors 
(see, e.g., [53155]). 

Fig. [9] illustrates the relative weight of the quasi-elastic and A production 
channels in the total inclusive 16 0(u e , e~) cross section. 



2 Notice that, compared to Eq. ([7j), the parametrization of the hadronic tensor in 
|52| involves an additional factor —1/2 in the coefficient of W$. 



1(3 




Fig. 9. Comparison between the CC quasi-elastic (solid line) and A production 
(dashed line) contributions to the total cross section of the process ie O(v e ,e~). 



Above E v ~ 0.5 GeV, the A production contribution becomes important; 
in particular, the ratio of the areas under the curves is of the order of 40 
%, thus showing the importance of including such a contribution in realistic 
estimates of the cross sections. To gauge the impact of A production on the 
measurements at /?-Beams, we have computed the averaged cross section 



(1(7)) = r , F ^ / 771 n t = QE,A 18 

for the fluxes reported in Ref . [20] . and constructed the ratio (a a) / (o~qe) as a 
function of the boost factor 7. The results are shown in Fig. [TUJ 

It clearly appears that the ratio of (cr)'s becomes larger and larger as the 7 
factor increases, until it becomes quite stable at 7 > 600 reaching the values 
0.52 and 0.37 for neutrino and antineutrino, respectively. This implies that 
the A production contribution to the total cross section is about 1/2 or 1/3 of 
the main contribution from CC quasi-elastic interactions. The presence of the 
plateaux is due to the fact that at higher energies the cross sections exhibit al- 
most the same energy dependence. Comparing with Fig. [T] we also understand 
the reason why the A production process contributes more to neutrino than 
to antineutrino cross section. For fixed 7, the antineutrino spectra is peaked 
at lower energies than the neutrino one, in a region where the contribution of 
the A is less important. In Appendix Owe list separately the values of (ctq E ) 
and (<ta) for the spectra of Fig. [TJ 
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Fig. 10. Ratio (ca) / (o~qe) as a function of the boost factor 7 at [3-Beams. The solid 
(dashed) line refers to neutrino (antineutrino). 

4.2 Contribution of higher resonances 

Beside the contribution of the A to the inelastic cross section, other reso- 
nances may be important, particularly for neutrino energies larger than ~ 2 
GeV. This so-called second resonance region contains three isospin 1/2 states: 
Pn(1440), -Di 3 (1520), 511(1535). New data on electroproduction have enabled 
the authors of [56] to estimate the vector form factors of these resonances. On 
the other hand, estimates for the axial couplings have been obtained from the 
decay rates of the resonances, exploiting the PCAC hypothesis. The strategy 
for computing their contributions to the cross section is the same used for A 
production, and we refer to [56] for a compilation of the form factors. Our 
results are shown in Fig. [TTJ 

As expected, even at the end of the energy spectra of Fig. [I] the contribution of 
the second resonance turns out to be marginal. For example, at E v = 1.5 GeV, 
the ratio of the cross sections is a a ■ cs u : o~d 1z '■ cp n = 1:0.12: 0.06 : 0.02. 



4-3 Comparison with other calculation and experimental data 

As for the quasi-elastic case, many papers have been devoted to the subject 
of resonance production. However, most of them focus on exclusive channels 
(a subject only marginally addressed in the following) and consider targets 
other than oxygen or carbon [571158] . For the sake of completeness, in Fig. [12] 
we compare our calculation for A production in 16 0(z/ e ,e~) interactions to 
the results of Ref . [4^] . The agreement between the two calculations appears 
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Fig. 11. CC resonance production via 16 0(^ e ,e~) interactions. The solid line shows 
the results for the A production, whereas the long dashed, dashed and dotted lines 
correspond to -D13, S\\ and Pu, respectively. 

to be better than in the case of quasi-elastic scattering. This feature may be 
ascribed to the fact that the approach of Ref. [49J includes the effect of the A 
propagation through the nuclear medium, neglected in our calculation. 

A production 

6 1 ■ 1 ■ . ■ 1 




0.0 0.5 1.0 1.5 



E v (GeV) 



Fig. 12. Comparison of our cross sections (solid lines) for 16 0(z^ e ,e ) and 
16 0(^ e ,e + ) to those of Ref. \49$ (dashed lines). 

To end this section, we compare our calculations to the preliminary results of 
the MiniBooNE collaboration for the cross section of the CC process ^ M + 12 C —■ 
fi~ + 7r + + X [59]. We should mention that a complete treatment of exclusive 
channels has not been carried out yet within our formalism. In particular, we 
are able to estimate incoherent 7r + production (in which the target nucleus can 
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be excited and/or broken up, leading to the production of pions) whereas we 
have not taken into account the production channel v + A — > fi~ + A + ir + , in 
which the target nucleus is left in the ground state (coherent pion production). 

As the CC one ir + production cross section is extracted from the measured 
ratio 



&CCQE 

using a specific model for CC quasi-elastic cross section, we prefer to calculate 
R from our model and then compare directly with Fig. 5 of [59]. We also 
include the contribution of the second resonance region. 

The following isospin decomposition can be applied to the pion production 
amplitudes A: [58] 



A(u e + p -y r + p + tt + ) = As, 

1 2 /2 

Aiyt + n^r + n + n + ) = -A 3 + ^-A u (20) 

fey 2 

A{u t + n -> r + p + 7i°) = — ^ 3 + 3 Ai, , 

where .A3 is the amplitude for the isospin 3/2 state of the irN system, pre- 
dominantly the A, and A\ is the amplitude for the isospin 1/2 state. From 
the above decomposition, we obtain the 7r + production cross section 



o{yt + p tt + ) = — o- A ++ + - (61 <r Pll + b 2 a Dl3 + 63 a Sll ) , (21) 

where the coefficients hi are branching ratios for n + production of the reso- 
nances (which we assume to be equal to 0.7 for Pn, Di 3 and Sn). 

The results of the calculation are shown in Fig. fflBl in which the experimental 
points are taken from [5H]; the dashed line refers to the A contribution only, 
whereas the solid line is obtained using Eq. 02 ip . 

The figure shows a fair agreement between experiment and theory. However, 
at energies above E v ~ 1 GeV our curves are below the data. This may be 
ascribed to non-resonant pion production [60J which has not been taken into 
account in the calculation. Another ~ 10% contribution could also come from 
coherent pion production [59]. A fully quantitative understanding of these 
effects deserves further studies. 
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Fig. 13. Comparison between the preliminary MiniBooNE data for the ratio 
R = a n + /acCQE in the process + 12 C — > yT + 7r + + X and the results of our 
calculation. The dashed line refers to the A contribution only, whereas the solid line 
is obtained considering also the contribution from the second resonance region (see 
Eq. m»- 

5 Conclusions 

In this paper we have investigated neutrino interactions with nuclei at energies 
up to ~ 1 GeV. In this energy region, the main contributions to the total 
inclusive cross section come from CC quasi-elastic and resonance production 
processes. To estimate the cross sections, we employed a formalism based on 
IA, applicable to both quasi-elastic and inelastic processes. 

For the quasi-elastic 16 0(V e ,e~) interactions, we compare various theoretical 
models, based on different descriptions of nuclear dynamics. We find that the 
simple RFGM overestimates the cross sections, with respect to the results ob- 
tained using a realistic spectral function to describe the strongly correlated 
bound nucleons. While dynamical FSI play no role in the total inclusive cross 
section, statistical correlations, arising from the effect of Pauli principle, fur- 
ther decrease the cross section, as the phase space available to the knocked out 
nucleons is reduced. Therefore, our best estimate of the inclusive total cross 
section is obtained using the LDA spectral function and including Pauli block- 
ing. These results are in fair agreement with those discussed in the existing 
literature. 

Although the formalism we have employed can be readily applied to a variety 
of different targets, for which realistic spectral functions are available, most 
of the available quasi-elastic neutrino data have been taken at low energies, 
where IA is not expected to provide accurate predictions. Higher energy data 
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will soon be provided by K2K SciBar experiment [61] and Minerz/a |17j . 

As for the inelastic cross sections, we focused on resonance production, in- 
cluding the contribution of both the A and the second resonance region. The 
calculated inclusive resonance production cross sections show that, as in the 
elastic case, the simple RFGM appears to be unable to accurately describe 
these processes. As expected, the CC resonanace production to quasi-elastic 
cross section ratio becomes increasingly important at E v > 0.5 GeV. At ener- 
gies below E u ~2 GeV, the contribution of Pn, -D 13 and Sn to the total cross 
section turns out to be small, being ~ 10%. 

Even if we do not consider exclusive channels explicitely, we provide an es- 
timate of the one ir + production through resonance excitation. Comparison 
between the calculated ratio R = joccQE to the preliminary results from 
the MiniBooNE collaboration shows a fair agreement between theory and data 
at E v < 1 GeV. On the other hand, at higher energies, a discrepancy emerges 
that may be removed once the contributions of non-resonant and coherent 
pion production are properly added to our results. 
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A Integration limits 

Let us analyze in detail the integration limits involved in the calculation of 
the quasi-elastic and anelasic cross sections. The integration in momentum in 
Eq. ([SD can be cast in the form 

d 3 p = p 2 dp d cos 7 d(j) , (A.l) 

where 7 is the angle between p and q and the integral on the azimuthal angle 
4> can be readily done, yielding a factor 2 w. Denoting by W the invariant mass 
of the final hadronic state produced in the interaction, cos 7 can be written in 
the form: 
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cos 7 = 



-w 2 - 


P 


2 _ 




q 


2 


+ (v + E p ) 2 


2 


P 




q 





or, equivalently, 

s + M 2 _ 1 -W 2 -2(v + M A )E A _ 1 n . 

cos 7 = o i i i i ' A - 2 

2 |p| |q| 

where M A _i = M A — + E, s is the squared center of mass energy and 
E\_ x = M\_ x + |p| 2 . The quantity cos 7 has to satisfy the constraints —1 < 
cos 7 < 1, leading to the lower and upper bounds for |p|: 



iPl = 1T 



1 

2s 



A|q| ± 



Mi 



A- 



AsM\_ 



1/2 



(A.3) 



with A = s + M\_ x — W 2 . The upper limit of the E integration can be found 
requiring the argument of the square root entering the definition of Ip^ to be 
non-negative. This leads to A > 2 \fs M A -i and finally to the bound 



En 



(W - m N ). 



(A.4) 



A.l Quasi-elastic case 



Up to now the integration limits we have found are completely general. The 
quasi-elastic case is recovered once we put W = m^. Equation (jHJ) can be cast 
in a simpler form once we use the energy conserving 5-function to perform 
the integration on cos 7 in Eq. (1A.1I) . We have to evaluate the Jacobian of the 
transformation 

d , 



m 2 N ) = 2 |p| |q| 



d cos 7 

leading to the final formula 



|^ = ^ [d P dE\ P \P(p,E)^. (A.5) 
dVldE e |q| J F IHI VH ' ' d£ldE e v ' 



A. 2 Resonance production 



For the case of resonance production, the variable W is the invariant mass of 
the resonance. It is more convenient to express Eq. (lA.ip in terms of |p| and 



23 



W instead of |p| and 00S7. It can be easily done with the help of Eq. (1A.2j) . 
implying 



2tt 

d 3 p = p 2 dp d cos j d(p =>- — 7 W dW Ipl dp . (A. 6) 

|q| 

The cross section for resonance production can then be written in the following 
form: 



d a IA 111 f f , f d (Jclcm 



dQdEf 



W th v- 



where d 2 a e \ era /dVldEi has the same form as in Eq. ([6]) but the hadronic tensor 
involves the Breit-Wigner factor of Eq. ( fl6l) and the structure functions are 
those specified in Ref. [52] . 



B Kinematical factors 



Here we report the kinematical factors appearing in the general formula for 
the cross sections, after tensor contraction (Eq. ([6])): 



2 

^2 



A 1 


= m l 


{k-k' 
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A 2 


= (k- 


p) (k' ■ 


p) 


A 3 


= (k- 


p) (k' ■ 


9) 


A, 


= (k- 


q) (k' ■ 


q) 


A, 


= (k- 


p) (k' ■ 


9) 



? A, 



2 m% 



2 ' 

N 



where 



(B.l) 



k=(E u ,k) k' = {E e ,k') 

P=(E\p\,p) Q= (p t q) . 

Notice that, due to the simmetry of the problem, some of the scalar products 
considerably simplify. In fact, using a reference system in which the vector q 
is along the z-axis and scattering takes place in the x^-plane, one can express 
the scalar products in terms of momentum component in that plane and, using 
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the symmetry of the problem, neglect terms linear in p x . The evaluation of 
the remaining contributions is straightforward. 

In the limit in which q —> q and p —> p, the coefficients A 2 -5 reduce to kinemat- 
ical factors of free elastic interactions. In particular, depend quadratically 
on the lepton mass, thus vanishing in the limit of massless leptons. This is no 
longer true in our formalism, in which, even in the limit of massless leptons, 
^4,5 are different from zero. The relative contribution of each term appearing 
in Eq. (jSJ) is shown in Fig. IB. II in which we plot the quasi-elastic cross section 
for electron neutrinos on l6 0. 




E v (GeV) 

Fig. B.l. Contributions of the different structure functions, defined in Eq. iB. 
to the quasi-elastic process o~( ie O(v e , e~)X). The values given in the shaded box 
correspond to the areas under the different curves. 

While the main contributions come from Ax_ 3 , non-vanishing A 45 give an 
additional 4% to the cross section . For heavier lepton masses, the contribution 
of A^z, is increasingly larger, as expected for free nucleon interactions. 



C Form factors for quasi-elastic scattering 

The hadronic current appearing in the definition of the hadronic tensor in 
Eq. ([3]) has the V — A structure J M = Jy — J A . The most geenral expressions 
of the vector and axial contributions read 

Jy = Jy + = U(p')T v T + U(p) 

J» = J^ = U(p>)T» A T + U(p) , 

where p (p') are the momenta of the incoming (outgoing) nucleon and T v and 
are given by: 
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2m N 



q"F s 



1 A 



In the definition of Jy and J a, r + is the isospin raising operator acting on the 
isodoublet U = (n,p) T and F\, F 2 , Fg, Fa, Fp, Ft are form factors depending 
on the squared four- momentum transfer q 2 . The CVC hypothesis allows one 
to relate the weak form factors to the eletromagnetic form factors and to 
constrain F$(q 2 ) = 0. On the other hand, the PCAC hypothesis relates Fp(q 2 ) 
and FA^q 2 ) in such a way that Fp(q 2 ) = 2m^ F^(g 2 )/(m^ — q 2 ). Finally, as a 
second class current is incompatible with the Standard Model of electroweak 
interactions [62J, we set Fx(q 2 ) = 0. 

To express the structure functions in terms of form factors we need to add 
the VV, AV and A A contributions to the hadronic tensor and equate the 
resulting sum to the decomposition (J7J). Thus, we obtain: 



Wi = 2 



~ {F 1 + F 2 f + (2m 2 N - q ^\F 2 A 



Fl 



4m 



W 2 = A 
W 3 = -4 (F 1 + F 2 ) F A 



,F 2 2 + F 2 

N J 



w 5 . 



F 1 F 2 + \2m 2 N + ' 



2 



F 2 



Defining the electric (Ge) and magnetic (Gm) form factors 
1 



G f 



i - -K 



Gm — 4.71 Ge 



(CI) 



and r = q 2 /Am 2 Nl we finally get 



F! 

F A 



1-t 



(Ge-tG 



M , 



1.26 



1 - r 
1.28 



i - -K 



1 _ j^ 2 V1.27 

0.14 



-Ge + Gm) 
Fa 



26 



D Averaged cross sections 



In Table (ID. II) we show the results of a for the quasi-elastic and A production 
induced by neutrino and antineutrino, for the two /3-beams scenarios described 
in the paper. 





7 


a QE (lO" 38 cm 2 ) 


cta (10~ 38 cm 2 ) 




100 


3.3 


0.28 


v e 


60 


0.34 


0.002 




250 


6.55 


2.9 




150 


1.35 


0.21 



Table D.l 

Average cross section a (see Eg. U8\) ) for quasi-elastic scattering and A production 
induced by neutrino and antineutrino. The two different sets of 7 values correspond 
to the fluxes shown in Fig. d 
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